Decoherence as Decay of the Loschmidt Echo in a Lorentz Gas. 
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Classical chaotic dynamics is characterized by the exponential sensitivity to initial conditions. 
Quantum mechanics, however, does not show this feature. We consider instead the sensitivity of 
quantum evolution to perturbations in the Hamiltonian. This is observed as an atenuation of the 
Loschmidt Echo, M(t), i.e. the amount of the original state (wave packet of width a) which is recov- 
ered after a time reversed evolution, in presence of a classically weak perturbation. By considering 
a Lorentz gas of size L, which for large L is a model for an unbounded classically chaotic system, 
we find numerical evidence that, if the perturbation is within a certain range, M(t) decays expo- 
nentially with a rate 1/t$ determined by the Lyapunov exponent A of the corresponding classical 
dynamics. This exponential decay extends much beyond the Eherenfest time tE and saturates at a 
time t s ~ A -1 ln(jV), where N ~ (L/a) 2 is the effective dimensionality of the Hilbert space. Since 
quantifies the increasing uncontrollability of the quantum phase (decoherence) its characterization 
and control has fundamental interest. 



Chaos justifies the observed macroscopic irreversibil- 
ity within the reversible laws of Classical Mechanics. One 
of its characteristic features is the exponential divergence 
of trajectories corresponding to nearby initial conditions, 
which leads to deterministic unpredictability. However, 
quantum dynamics exhibits insensitivity to initial con- 
ditions [|l| and hence prevents a dynamical definition 
of quantum chaos. Therefore, quantum signatures of 
chaos in systems with a chaotic classical equivalent are 
searched in their steady-state properties Q such as spec- 
tral rigidity || , wave function morphologies Q , and de- 
caying parametric correlations of observables [5[. Early 
attempts Q to address unitary quantum dynamics did 
not clarify the connections with the dynamical classical 
concepts of chaos, but the inclusion of interactions with 
a dissipative environment was expected @ to produce an 
entropy growth controlled by chaos. In a purely Hamil- 
tonian problem, quantum reversibility can be monitored 
through the amount of revival of a local density excita- 
tion, upon time reversal of its unitary quantum evolution, 
i.e. the Loschmidt Echo By considering a Lorentz 
gas where the reversed evolution is disturbed by a static 
perturbation, we find a Loschmidt Echo that attenuates 
exponentially with a rate associated to the chaos of the 
classical system. Under appropriate conditions, the dy- 
namical instability of the classical system translates into 
quantum phase unpredictability (decoherence) and the 
classical Lyapunov exponent becomes the quantum irre- 
versibility rate. 

A hypersensitivity of time reversal to perturbations 
was observed in recent NMR experiments on many-body 



spin systems [|| [|l0|. In essence, these experiments Q 
p"2| involve the creation of a local density excitation rep- 
resented by a state |0) which evolves under a many-spin 
Hamiltonian H. The Loschmidt Echo (LE) is the proba- 
bility to return to the initial state when a Hamiltonian 
evolution for a time t is followed by an identical period 
of imperfect reversal of that evolution, achieved by the 
transformation 7i — » — (7i + £), i.e.: 



M(t) 



| exp[i(H + T,)t/h] exp[-iHt/h] |0) 



(1) 



where £ is a constant (or quasi-static) Hcrmitian pertur- 
bation representing the imperfection in the Hamiltonian 
reversal. Notice that in these experiments — ln[M(£)] is 
a measure |^ of the growing entropy. A striking finding 
|]l0f was that, for small E, the decay of M(t) becomes in- 
dependent of E, being proportional only to the dynamical 
scales of TC. The hints that chaotic systems may be un- 
stable under fluctuating perturbations [Q and that many- 
body systems could be assumed to be intrinsically chaotic 
|fl3f suggested us the hypothesis that the most relevant 
parameters of the Hamiltonian dynamics, the Lyapunov 
exponents A, might control the decay of M(t). More- 
over, a semiclassical analysis, assuming |l4) a one-body 
classically chaotic Hamiltonian with a perturbative po- 
tential representing a quenched disorder predicted a LE 
decaying exponentially with a rate 1/t^ = A. These pre- 
dictions, based on various approximations, encouraged us 
to perform exhaustive numerical experiments. 

We consider a Lorentz gas, i.e.: a particle in a square 
billiard of area L 2 where we fix an irregular array of ./V 
circular scatterers (impurities) of radius R. A particular 
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FIG. 1. Wave packet evolution in a system with L/a — 200, 
i?/a = 20 and c = 8iv x 10~ 2 . a) Initial state at t = 0. 
The velocity points to the left, b) State evolved with the 
unperturbed Hamiltonian for a time t = 30H/V. c) State 
evolved with the perturbed Hamiltonian (a = 0.04) for a time 
t = 30h/V. d) State evolved from that depicted in panel b) 
with the perturbed Hamiltonian for a time t — —30H/V. The 
square of the overlap (Loschmidt Echo) between the states of 
panels a) and d) is M(t) = 0.09, same as that between the 
states of panels b) and c). 

realization of such system is represented in Fig. |]-a where 
the length, in a minimal unit a, is L = 200a. The scat- 
terers concentration is c = NttR 2 /L 2 = 8n x 10~ 2 . The 
minimum distance allowed between two scatterer centers 
is 3i?. This requirement, together with the imposed pe- 
riodic boundary conditions, prevents the geometrical lo- 
calization. 

The Lyapunov exponent of a Lorentz gas should be: 
A = (3v/£, where I is the collision mean free path, v the 
particle velocity and /3 ~ 1 a geometrical factor [|l5| de- 
pending logarithmically on £. In our system, we estimate 
£ ~ L 2 /(2NR) - -kR/2. The computation of the distribu- 
tion of distances between collisions gives a shifted Poisson 
distribution whose mean value is I. The Lyapunov ex- 
ponent is obtained |]l6f from the average logarithm of 
the distance dt, after a time t, between two classical tra- 
jectories initially separated by a small distance do- The 
condition of smallness for do is that dt <C R. The longer 
the time the more precise the estimation of A. By ne- 
glecting correlations in the position of the impurities, we 
obtain the estimation: 







ds 



exp(— s/£) 
eexp(-R/l) 



The outer integral accounts for the distribution of free 
paths, s, and the internal one is the distribution of impact 
parameters x. Numerical X£/v verify Eq. (Q). For the 
present case of c = 87r x 10 -2 , i.e. I — 96a, we got (3 
= (2.1 ±0.1). 

The perturbation parameter a controls the distortion 
of the diagonal components of the mass tensor , m XtX = 



m (l + a), and m y 



i /(l + a), with m the isotropic 



unperturbed mass. This perturbation is inspired in the 
effects of a slight rotation of the sample in the related 
problem of dipolar spin dynamics flT7| ], which modifies 
the mass of the spin wave excitations. It is written as 



S(a) 



1 + a 2m n 



Py 
t— — 

2m„ 



(3) 



dx ln[l 



RVl 



(2) 



For a fixed initial position and velocity we find nu- 
merically that two evolutions with slightly different mass 
tensors lead to trajectories whose distance in phase space 
grows exponentially with the same Lyapunov exponent 
that amplifies initial distances: i.e. the classical dynam- 
ics of this system is equally sensitive to changes in the 
Hamiltonian as to changes in the initial conditions [Q. 

The Hamiltonian operator is obtained by a lattice dis- 
cretization on a small scale a. The evolution is calcu- 
lated using the Trotter-Suzuki (TS) algorithm [jl9|. Its 
basic idea is that, by choosing a short time unit r, the 
evolution can be approximated as the finite product of 
evolution operators where each is solved analytically. In 
the lowest order, 



U(t) = exp[mT/h} ~ U{t) = JJcxp[iW fc T/fi], (4) 

k 

where Ti- — Y^k ^ fc - The conceptual virtue of this method 
is that U(t) is always unitary. For an approximation of 
order n, the difference between U(t) and U(t) is order 
t" +1 . Our choice of n = 4 and t = O.lh/V = 0.2m o a 2 /h 
allowed us to obtain high accuracy even for times an order 
of magnitude larger than needed. 

Let's consider a typical system with R/a = 20, c = 
8tt x 10~ 2 and size L/a = 200 . The initial state is a 
Gaussian wave-packet of width a = 3a and wave number 
k x = — 0.7 /a, k y = shown in Fig. [p-a. All further simu- 
lations will use k x = k y = 3ir/8a. Since kR S> 1, compar- 
ison with semiclassical calculations is justified. We took 
27r/fc <C £ = (96 ± 2)a, in order to minimize Anderson's 
localization effects. The density resulting from a typical 
evolution for a time t = 30^/1^ |^0| is shown in panel b). 
An evolution with the perturbed Hamiltonian (a — 0.04) 
for the same time is shown in Fig. ^-c. In both pan- 
els the classical trajectories corresponding to three ini- 
tial positions are shown for reference. While densities 
look similar to the eye, M(t) is about 0.09 indicating the 
relevant role of the quantum phase in the attenuation of 
M(t). Panel d) shows the LE formed by the perturbed 



2 



0.1 



0.01 



1E-3 



1E-4 



1E-5 



\ \\ 1 £. - 




\ ? 8 




: 1 4 




u\ n 




•A 0.0 0.7 1.4 2.1 

\\ 1/L 2 [10V] 


i Jo 


1 I-r U^-, , r-^C 



t E 150 



300 



450 600 

t [fiN] 



750 900 



FIG. 2. In bold line, M(t) for an individual system with 
L/a — 800 and a = 0.024. The Ehrenfest time is shown with 
an arrow. In dashed line, shown for reference, exponential 
decay with the calculated Lyapunov exponent of the classical 
equivalent. In solid lines, average M(t) for different system 
sizes L/a =200, 400 and 800 and the same perturbation. In 
the inset M x is shown for L/a = 200, 300, 400, 500, 600 
and 800 as a function of (a/L) 2 . The straight line is the best 
fitting, with ¥„ = (5 ± 0.2)(a/L) 2 



backward evolution of the state in panel b). Analyzing 
M(t) in different realizations we observe that, after the 
initial transient, M(t) fluctuates around an exponential 
decay, with a characteristic time r$. Notably this expo- 
nential decay persists much beyond the Ehrenfest time 
t E ~ \ \n[2R/\ F ] ~ 40 h/V, which is fixed by the local 
scale of the potential fluctuation 0. In a Lorentz gas, 
in contrast with the usual case of chaotic cavities, t E is 
independent of the system size Finally, M(t) sat- 

urates in a time scale related to system size. A typical 
curve for L/a = 800 is shown with a bold line in Fig.(||). 

In principle, any given precision of can then be ob- 
tained studying a system large enough, however this eas- 
ily becomes computationally expensive. Alternatively, 
one can resort to the ensemble averaging of the observ- 
able M(t), which reduces noise and defines with the 
same precision in much smaller systems. This is exem- 
plified in Fig. g, where we present the ensemble average 
M(t) for billiards of three different sizes of L/a =200, 
400 and 800 and fixed a = 0.024 and c = 8tt x 10~ 2 
which show the same exponential behavior with a pro- 
gressively expanded range. Similar plots are obtained for 
the other parameters with an exponential M(t), showing 
that the ensemble average does not introduce any spuri- 
ous effect in the decay. Previous attempts to character- 
ize quantum chaos flq] considered the limiting value 



for the Loschmidt Echo at asymptotically large times. 
It was proposed oc 1/N, where iV is the num- 

ber of energy levels appreciably represented in the initial 
state. In our case N ~ (L/a) 2 . The numerical results 
verify a complete independence of on the pertur- 
bation. The calculated for various sizes are shown 
in the inset of Fig. (g) for a fixed a = 0.024 . The 
predicted relationship is shown with a straight line ob- 
taining Mqo = (5-0 ± 0.1)(a/L) 2 . For all the sizes we 
kept c (and hence £, A and t E ) fixed and we verified that 
T0 does not depend on L/a. These results imply that 
the Lyapunov exponential behavior persists up to a time 

A representative dependence of on the parameter a 
is obtained considering the smallest sample size compat- 
ible with a good observation of the exponential, in this 
case L/a = 200. Fig. || shows M(t) averaged over 100 re- 
alizations that contain at least a scatterer in the classical 
trajectory of the wave packet. For a > a c ~ 0.02, all the 
exponential decays coincide, within the numerical error, 
with the Lyapunov decay associated with the classical 
system, shown with a thick line for comparison. The ini- 
tial perturbative parabolic decay, M(t) ~ 1 — b (at) 2 with 
b ~ 0.37 (V/h) , has a characteristic scale which does 
depend quadratically on the perturbation strength and 
prevents the curves to lie superimposed. We show in the 
inset of Fig. g the numerical values of t^ 1 extracted from 
the exponential part. It shows an initial quadratic depen- 
dence on a and a cross-over at cv c ~ 0.02 to a saturation 
at a value close to the classical Lyapunov exponent. 

The semiclassical analysis predicts a universal behavior 
provided that the perturbation is strong enough to be 
quantically significant, but weak enough to be classically 
irrelevant. The first condition implies that the length 
£ — vt required for an important change in the phase 
as consequence of the perturbation must be shorter than 
the distance associated to the Lyapunov exponent v/X. 
A rough estimation for our perturbation Sm x . x — am Q 
is: 



5E~ k 



irm a 



ah(l — cos[fcaj) 



< 1/A 



(5) 



which would predict a critical value a c ~ 0.13. However, 
we see the universal exponential even for perturbations 
a > 0.02 implying a critical value an order of magnitude 
smaller than Eq. (||) . This suggests that either this esti- 
mation or the condition required in Ref. Jli] is too strong 

II- 

The condition for a classically weak perturbation 
means that it must not modify the system's global prop- 
erties. Otherwise, even in absence of chaos (no colli- 
sions) the perturbation would spread out the classical 
trajectories. In our case, this M(t) fits a Gaussian de- 
cay. If at the time 1/A, required for the chaos to set 
in, the overlap has already decayed by cause of the per- 
turbation, the exponential decay would not be observed. 
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FIG. 3. Average M(t) for a system with L = 200a, 
R = 20a, which makes £ — (96 ± 2)a. The values of a are, 
from top to bottom: 0.006, 0.009, 0.012, 0.018, 0.024, 0.03, 
0.04, 0.05 and 0.06. The thick line corresponds to an expo- 
nential with the calculated numerical Lyapunov exponent of 
the system. For long times M(t) saturates to a finite value, 
Moo. In the inset the inverse characteristic decay times are 
shown, evidencing the regime where the decay of M is given 
by the intrinsic properties of the system. The dotted is the fit, 
tT 1 — (130 ± 6)a 2 V/h, while the dashed line is the classical 
Lyapunov exponent. 



suits through their effect on the wave function's phase. 
In this sense, the definition of classical chaos in terms 
of hypersensitivity to initial conditions and perturba- 
tions translates, in the quantum world, as a sensitivity 
of phases to perturbations in the Hamiltonian, i.e. deco- 
herence. Since we restricted ourselves to a Hamiltonian 
problem, its solution clarifies how chaos limits our con- 
trol of dynamics even within the reversible framework of 
Quantum Mechanics of closed systems. This concept has 
deep consequences for the diverse fields where decoher- 
ence must be minimized, such as classical waves prop- 
agation p3[ , mesoscopic and molecular electronics |24|] 
and quantum information [ p5| . In view of the entropic 
meaning of — ln[A/(<)] the interest for foundations of Sta- 
tistical Physics Q is clear. 

In summary, we have presented numerical evidence 
supporting the fact that decoherence, as measured by the 
exponential decay of the Loschmidt Echo, is controlled by 
the same parameters that govern classical chaos. This 
renders M(t), defined in Eq. [|, as a very practical tool 
to study dynamical quantum chaos. 

We acknowledge the Argentinian Supercomputer Cen- 
tre and the UNC Department of Computer Science. 
DAW and FMC receive support from CONICET and 
SeCyT-UNC respectively. HMP and PRL are affiliated 
to CONICET. This work received financial support from 
CONICET, ANPCyT, SeCyT-UNC, Fundacion Antor- 
chas and the french-argentinian ECOS-Sud program. 

* Correspondent author. E-mail: 
horacio@famaf.unc.edu.ar 



This sets an upper bound for the perturbation from 
M(t = 1/A) = cxp[-6(a/A) 2 ] < M^. Replacing with 
the parameters of the system shown in Fig. (^) one gets 
the condition a < 0.3, consistent with our results. 

Therefore, we have shown that in a wide range of pa- 
rameters M(t) in any large enough individual system de- 
cays exponentially until an asymptotic value is reached. 
The characteristic decay time does not depend on the 
perturbation, but rather on the intrinsic properties of 
the Hamiltonian. The range of perturbation parameters 
with an exponential decay results to be much broader 
than hinted by the previous semiclassical analysis. From 
our numerical results one can infer that infinite systems 
not showing Anderson's localization should present an 
unbounded exponential decay. The Lorentz gas repre- 
sents a broad class of chaotic systems, where it is ex- 
pected the same overall behavior. It has the additional 
advantage over other chaotic cavities that here the sat- 
uration time (t s ~ ^ln[L/a]) can be made arbitrarily 
longer than the Ehrcnfcst time. This manifests the quan- 
tum nature of the observed effect. The simplicity of our 
perturbation evidences a very remarkable property: there 
is no need of chaoticity or stochasticity in the "environ- 
ment" to introduce irreversibility. Fig. [l] manifests that 
the irreversibility produced by small perturbations re- 
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